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Introduction 


As  demonstrated  in  Lowengrub  [ $ ] 


reduce  to  the  problem  of  determining  functions 


satisfying  Abel  type  integral  equations  of  the  type 


is  Holder  continuous 


In  this  paper  we  show  that  systems  of  the  type  (1.1)  and 


(1.2)  can  be  reduced  to  the  determination  of  a matrix  function 


satisfying  certain  growth  conditions  at 


is  a coefficient  matrix  with  elements  linear 
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combinations  of  ar,p,v  and  6.  The  system  (1.3)  is  a coupled 
system.  We  effect  a linear  uncoupling  of  this  system  by  intro- 
ducing certain  similarity  transformations.  The  matrices 
associated  with  these  transformations  are  explicitly  computed 
and  hence  exact  solutions  are  derived,  in  physical  applications, 
such  as  the  determination  of  the  stress  field  in  an  inhomogeneous 
body  containing  flaws,  the  relevant  physical  quantities  are 
expressed  in  terms  of  the  matrix  function  $(z).  One  need  not 
actually  solve  for  ^ and  cp2  in  (1.1)  and  (1.2).  The 
functions  §^(z)  and  §2  (z)  are  defined  (say  in  the  case  p = 1) 


•i‘»>  ■ ife  I 2^ 

a (t  - z)u 


Art!. 

where  R(z)  = [(z-a)(b-z)]  2 


These  functions  must  be  defined 


on  appropriate  branches.  Analogous  representations  are  intro- 
duced for  p * 1 . 

Section  2 of  the  paper  thoroughly  analyzes  the  case  p = 1 
while  in  section  3 we  choose  p = 2.  These  are  the  two  cases  of 
physical  interest.  In  section  4 we  consider  some  explicit 
examples:  (i)  ar(x)  =a,  p (x)  = p,  y(x)  = y and  6 (x)  =6  with 
\i  = H,  P = 1 and  of,p,v  and  6 constant;  (ii)  o(x)=p, 
p (x)  =5,  y(x)=a  and  6 (x)  = with  li  = *s,  P = 2,  and  a. 


yu  . 


P constant.  The  final  section  demonstrates  how  general  dual 
relations  (given  in  terms  of  the  Erdelyi -Sneddon  operators  of 
fractional  integration  - see  Erdelyi-Sneddon  [ l3)  may  be 
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2.  The  first  generalized  Abel  svatem. 

In  this  section  we  consider  the  generalized  Abel  system  of 
equations. 


pxqp1(t)dt  bqp2  (t)dt 

a(x)J  ~U-t)V  + e(x)  J = fl(x)'  X€  (a'b)  C 2 • 1 ) 


b cp-^tjdt  x cp  (t)dt 

Y(X)J  7t^3“  + 6(x)J  - f2(x)*  x € (a,b)  (2.2) 


x a 

where  a,$,Y,  and  6 satisfy  conditions  to  be  specified  later. 


However,  we  do  assume  that  f^  and  f^  are  Holder  continuous 


on  (a,b). 

As  in  Sakalyuk  [(?],  we  define  the  sectionally  analytic 
functions 


$i(z)  = r7z)J  Tt-z)»  dt'  i = l,2. 


(2.3) 


1^ 


where  R(z)  = [ (z-a)  (b-z) 3 z and  the  function  is  defined  by 


some 


branch.  If  qp*(t)  satisfies  q,  (t)=  q>  (t)  [ (t-a)  (b-t) 
where  e>0  and  cp?(t)  is  Hblder  on  [a,b],  then  ^(z)  is 


analytic  in  the  plane  cut  along  [a,b].  Moreover,  the  boundary 


limits  (x) , where 


*±(x)  - lim  *.  (z)  a < x < b 

* Z"*X  J- 

Im (z)>0 


and 


* . (x)  = lim  i ( z ) , 


z-x 


a < x < b 


lm(z)<0 


1 


- 


6 


are  continuous.  In  addition. 


u-.l 

(z)  = o (|  z - a|^~  ) , 


^(z)  = 0 ( | z - b|  2 ) / z-*b 


(2.4) 


$.(z)  = Odzl-1)  , 


A simple  calculation  verifies  that 


$*(x) 


_i_  r rb  cpi(t)dt  -i 

R(x)le  J (x  - t)H  + J (t  -x)^  J 


(2.5) 


x € (a,b) 


*.(x) 


1 rpXcPi(t)dt 

‘ R(x)U  (x-  t)^  + 6 


. b m.  (t)dt  _ 
MTTif  i 

J„  (t  - x)  ^ J 


It  immediately  follows  that 


xep.(t)dt  e^nif^(x)  +«"(x) 

I (x-  t)M  = [ e^TTi_  1 JR  (x)  / 


X € (a,b) 


(2.6) 


(2.7) 


pb  ?i(t)dt  p **(x)  +eMTTi*“(x)  ^ 

j (t-x)U  = - 1 e^-1 ><x)  * 


(2.8) 


i - 1.2  . 


Substitution  of  (2.7)  and  (2.8)  into  (2.1)  and  (2.2)  yields 


the  boundary  condition  on  a < x < b , 


a (x)eUrri$*  (x)  - 3 (x)  $2  (x)  + a (x)  (x)  - 3 (xJe^S”  (x)  = F1  (x)  (2.9) 

-\(x)$^(x)+6  (x)eurri$2  (x)-v  (x)eUTT1§^  (x)+6  (x)  $ ~ (x)  = F2  (x)  (2.10) 


where 


Fi(x)  = fi(x)  (e2tiTT1  - 1)/R(x) 


This  substitution  then  reduces  our  problem  to  determining  two 
sectionally  analytic  functions  $1(z),  $2(z)  satisfying  the 
growth  conditions  (2.4);  that  is,  solving  a coupled  Riemann- 
Hilbert  boundary  value  problem.  Once  we  have  determined  (z), 
(i  = l,2),  the  functions  cp^  and  cp2  are  obtained  by  solving 
the  Abel  integral  equations  (2.7)  and  (2.8). 

For  convenience  we  introduce  the  following  matrix  notation; 


§ (z)  = (^(z)  ,$2(z))T  and  F (x)  = (F1  (x) , F2  (x) 

so  that  the  boundary  conditions  (2.9)  and  (2,10)  may  be  written 
in  the  form 


A(x)*  (x)  = - eMTT1A(x)  T(x)  +F(x)  , a < x < b (2.11) 


where 


A (x)  = (a_(x)),  i,  j = 1,2  with 


aH  (x>  = a (*)  ■ ” ' a12  (x)  = -p  (x) , a21  (x)  = -y  (x)  and 

a22(x)  = 6 (x)euni. 


(2.11a) 
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If  we  require  the  condition 
det  A (x)  = a (x)  6 (x)  e^771-  v (x)  g (x)  0,  a < x < b 


(2.12) 


then  the  matrix  A(x)  is  invertible  and  (2.11)  is  equivalent 
to  the  boundary  condition 


$+(x)  = -e^^G  (x)  $ (x)  +g(x) 


a < x < b 


(2.13) 


where 


G (x)  = A ^(x)A(x)  and  g(x)  =A  '*'(x)F(x). 


(2.14) 


It  is  necessary  for  us  to  determine  conditions  whereby  the 
coupled  Riemann-Hilbert  problem  can  be  uncoupled.  We  shall 
effect  a linear  uncoupling  of  the  system  (2.9)  and  (2.10)  by 
finding  a non-singular  matrix  P(z),  analytic  in  the  complex 
plane  (except  for  perhaps  a finite  number  of  poles)  with  a pole 
at  infinity  and  such  that  for  a < x < b 


P(x)G(x)P  (x)  = D(x)  = (d±j  (x) ) 


with  d12 (x)  = d21  (x)  = 0 . 

T 

Let  £(z)  = (E^ (z) * ^ (z) ) be  defined  as  follows: 


£(z)  = P (z) 9 (z) 


(2.15) 
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Note  that  £( z)  is  analytic  in  the  plane  (except  for  perhaps 
a finite  number  of  poles)  cut  along  a<x<b.  In  addition, 

£(z)  satisfies  appropriate  growth  conditions  at  infinity. 

Substitution  of  (2.14)  into  the  boundary  condition  (2.13) 
yields  the  uncoupled  set  of  conditions, 

£*  (x)  = -eUTTidii  (x)  £i  (x)  + k.(x)  , a<x<b  (2.16) 

i = 1, 2 and  k(x)  =P(x)g(x).  Thus,  this  procedure  reduces  our 
problem  to  the  determination  of  two  sectionally  analytic  functions 
^(z)  and  £ 2 (z)  satisfying  appropriate  growth  conditions  at 
infinity  and  the  boundary  conditions  (2.15).  The  solution  to 
these  Riemann  Hilbert  problems  is  well-known  (once  the  index  is 
determined)  (see  for  example,  Muskhelishvili  [£]). 

The  main  problem  is  to  explicitly  determine  the  matrix  P. 

We  first  compute  G(x)  from  (2.14)  . 

G (x)  = [det  A (x) ]_1  T (x) 
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where  T (x)  = [t^(x)]  with 

til  = t22  = o' (x)S  (x)  - g (x)  y(x) , t12  = -2ig  (x)6  (x)  sinun 

and  t21  = -2ia  (x)  v (x)  sin  utt  . 


A non-singular  matrix,  P(x),  for  which  P(x)T(x)P  ^(x) 
is  diagonal  exists  if  and  only  if  T(x)  has  two  linearly  inde- 
pendent eigenvectors.  In  this  case,  the  matrix  P(x)T(x)P  ^(x) 
has  the  eigenvalues  of  T(x)  as  its  diagonal  elements  and 
P (x)  has  as  its  rows  two  independent  eigenvectors  of  T(x). 

If  t^2  (x)  t^  (x)  ^ 0 then  the  eigenvalues  of  T(x)  are 
tn  (x)  + Jt  (X)t  (x)  so  that  for  the  matrix  P(x)  we  may 
take 


P(x)  = c(x) 


^(x)  - Vt12(x)t21(x) 


^(x)  Vt12(x)t21(x) 


(2.17) 


where  c(x)  is  any  scalar  function. 

If  t12 (x) t21 (x)  - 0,  then  T (x)  has  two  independent  eigen- 
vectors if  and  only  if  both  t^2  (x)  an<^  *"21  vanish. 

Since  T (x)  is  non-singular  this  can  occur  if  and  only  if 
either  6(x)=a(x)=0  and  y(x)@(x)^0  or  y(x)  =g(x)  =0  and 
& (x)a(x)  / 0.  In  this  case,  T(x)  is  just  a scalar  multiple  of 
the  identity.  For  most  applications,  t^2  (x)  t^  (x)  ^ 0 except 
perhaps  at  the  endpoints  x = a and  x = b.  Such  exceptional 
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cases  are  easily  handled.  In  order  to  simplify  further  work, 


we  assume  that  t^2  (x)  t^  (x)  ^ 0 . 

Let  P(x)  be  defined  by  (2.17)  so  that 


P(x)T(x)P  (x)  - 


tlx (x)  + Jt12  (x) t21 (x) 


tn(x)  - Vt12(x)t21(x) 


which  provides  the  desired  uncoupling  of  the  pair  of  coupled 
Riemann-Hilbert  boundary  value  problems.  The  scalar  c (x) 
is  chosen  so  that  P (x)  can  be  extended  to  a matrix  P(z) 
meromorphic  in  the  plane  with  a pole  at  infinity.  This  enables 
us  to  explicitly  determine  and  S2  and  hence  and  $2* 

Inversion  of  (2.7),  (2.8)  yields  our  original  unknown  functions 


cp.  and  cd, 

1.  % 4 


An  example  of  this  analysis  appears  in  section  4. 


3 . The  second  generalized  Abel  system. 


We  next  consider  the  system 


x cp  (t)  dt 


1 cp2  (fc)dt 


r 1 r 

a(x)J  2 2 |i  + P(X)J  ~~2  2 

0 (x  -t  )M  x (t  -x 


(t  -x“) 


■ = f (x)  , 0 < x < 1 (3.1) 

I* 


1 cp  (t) dt  x cp7  (t)  dt 

(x)  f \ -2-m  +Mx)  V'  o 'n  = f2(x)  ' 0<X<1.  (3.2) 

x (t1 2-xV  J0  (X  -t  )M 
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The  interval  (0,1)  is  chosen  rather  than  (a,b)  in  order  to 
simplify  the  analysis.  It  is  trivial  to  extend  to  (a,b). 

The  case  p = 1/2  in  (3.1)  and  (3.2)  has  been  considered  by 
Lowengrub  [3  ] . 

Analogous  to  the  method  used  in  section  2,  we  introduce 
the  sectionally  analytic  functions  ^(z)  and  §2  (z)  defined 

by. 


$ . (z)  = (z^_ 

l 


cpi  (t ) dt 


(*2-tV  ’ 


i = 1,2 


(3.3) 


If  <p.(t)  satisfies  cp  . (t)  = cp*  (t)  t^[t  (1-t)  ]|J+e  ^ where 

X 1 1 

e = 0 and  cp?(t)  is  Holder  continuous  on  [0,1],  then 
§^(z)  (i  = l,2)  is  analytic  in  the  plane  cut  along  [-1,1]  and 
satisfies  the  following  conditions, 


and 


= o(| 

as 

z -»  1 

(z)  = 0(| 

as 

N 

i 

i 

$i(z)  = 0(| 

z|_1) 

as 

Z -»  00  . 

(3.4) 


Moreover,  the  limiting  values  are  continuous  funotions 

|x|  <1  except  perhaps  for  x = 0. 

For  each  of  the  functions  (z-1)^1  **,  (z+l)^  ^ , (z  - t)  ^ 

(z+t)  ^ we  take  as  the  branch  cut  that  line  lying  along 


I 


for 


and 


o , f 1 tp.  (t)dt  . x cp.  (t)dt 

^W-iU-xV-15  +.Mnir  !i_. 

(f -xV 


0 (x2-t2)M 


„ , 1 cp.  (t)  dt  . * m.  (t)dt 

2iMs  P li ^ ^-UTiir  ^1.1  


$.  (x)  = i (1— x ) . , 

i J ..2  2.u 

X (t  -X  )K 


+ e 


0 (x2-t2)u) 


whereas  for  -1  < x < 0 


§*(x)  = -i(l-x2)U-ls 


i ^(tiat  unl|.’t|iP1(t)at  1 

J 2 2 u J , 2 2,|j  | 

| xj  (t  "x  > o (x  ~ t )»J 


2 p1  cPi(t)dt 

§i(x)  = i(l-x2)M  2'  r 


X 


<t2-*V 


+ e 


um 


'p  ,cpi(t)dt 


‘o  (x2-tV 


It  should  be  observed  that  for  -l<x<0. 


± I 

$i (x)  = - §i (-x)  . 


. _ o > x <P-  (t)dt 

(x)  + 5 (x)  = 2 sin  htt(1-x  )M_  2|  — ^ — 7-  , 


0 (x  -t  )1 


and 


.-MTT1.+ 


1 cp,  (t)dt 


(3.5) 


(3.6) 


(3.7) 


(3.9) 


(3.10 


L§^(x)  +eMTTi$i(x)  =-2  sin  utt  (1-x2)  f - ~2~~~2~  (3.11 

'X  (t  -x  )u 


for  0 < x < 1. 

Substitution  of  the  above  into  the  original  set  of  integral 


equations  (3.1)  and  3.2)  reduces  the  problem  to  the  following: 


determine  two  sectionally  analytic  functions  (z)  and 
$2(z)  satisfying  the  conditions  (3.4)  and  the  boundary  values, 

Ca(x)$^(^)-e(x)e~UTT1$2(x)]  + [a(x)^-g(x)eMni$‘(x)]  = Fx  (x)  (3.12) 

0 <x<  1 

[-V(x)e_MTTif^(x)+6  (x)§2(x)]  + [6  (x)  *"  (x) -V  (x)  eUTT1  (x)  ] = F2  (x)  (3.13 

0 < x < 1 


where 


Fi(x)  = 2 sin  4rr(l  - x2)u_^fi  (x)  . 


In  matrix  notation  the  system  (3 . 12) - (3 . 13)  becomes 

A(x)$+(x)  = - B(x)$  (x)  + F (x)  0<x<l  (3.14) 

where  $ (z)  = (^  (z) , §2  (z)  )T,  F (x)  = ^ (x) ,F2 (x) )T, 

A (x)  = (a^  j (x)  ) with  ai;L(x)  = a (x) , 

al2  (x)  = -p  (x)e"urri,  a21  (x)  = -Y(x)e-|aTT:L  and  a22  (x)  = 6 (x)  , 
and  B (x)  = A (x)  . (3.14a) 

In  order  to  determine  $(z),  boundary  conditions  for 
$ (x)  must  be  extended  to  all  of  [-1,1].  It  is  clear  from 
(3.9)  how  this  extension  is  to  be  performed. 

In  particular,  we  obtain  the  system 


A (x)  i+  (x)  = - B(x)$(x)  +F(x) 


-1 < x < 1 (3.15) 


where  A (x)  = 


B(x) 


A (x) 

0 <X<  1 

-A  (-x) 

-1  <X<  0 

B (x) 

0 <x<  1 

-B(-X) 

-1  <x  < 0 

(3.16) 


KXC 


15 


and 


A 

F(x)  = 


F(x) 


0 <x  < 1 


F(-x)  -1  < X < 0 


As  in  section  1,  we  assume  that  K (x)  is  invertible  for 
-1 < x < 1 . The  expressions  in  (3.16)  imply  that  it  suffices  to 
assume  that  A (x)  is  invertible  for  0<x<l,  or  equivalently, 
that 

a (x)6  (xje2^1-  y (x)  p (x)  ^ 0 0<x<l  . 


This  gives  us  the  system, 

§+(x)  = -G(x)$(x)  + 6 (x) 


-1<X<1  (3.17) 


G(x)  = A ^(xjBlx)  = A ^xjAfx) 


where 
and 

Thuft  we  must  now  find  a sectionally  analytic  matrix,  $ (z) 


g(x)  = A-1  (x) F (x)  . 


satisfying  appropriate  conditions  at  1 and  -1,  a growth  condition 
$ (z)  = 0 (| z| -1)  at  infinity,  and  the  boundary  values  (3.17). 

We  seek  a linear  uncoupling  of  (3.17);  that  is,  a non-singular 
matrix  #(x)  such  that  P (x) G (x)P-1 (x)  is  a diagonal  matrix 
for  -1  < x < 1. 

A 

We  first  compute  G(x)  on  0<x<l  and  obtain, 

G(x)  = [det  (A  (x) ) ]-1T  (x) 

where  T (x)  = C t^  (x)  ] with  t^  (x)  = t22  (x)  = 6 (x)o-  (x)-y(x)  p (x) , 
t12  (x)  = -2i  sin  pttS  (x)  p (x)  and  fc21  (x)  = -2i  sin  |jTTa  (x)  y(x)  . 

If  we  employ  the  arguments  of  section  1,  we  see  that. 


L- 


' ' ' 


WIPHUlPWilH  ' 


p (x)  = c (x) 


12(X) 

-Vti2 (x) t21 (x) A 

12<*) 

Vt12(x)t21(x)  j 

• 

(3.18) 

= + 2i  sin  (jn  c (x) 


-&{x)g(x)  —*Jot  (x)  g (x)  V (x)  6 (x) 
-6(x)p(x)  Ja (x) g (x) y(x)6 (x) 


(provided  t12(x)t21(x)  jt  0) 

produces  the  desired  uncoupling  of  the  boundary  conditions  (3.17) 
for  9 < x < 1. 

In  particular, 

An  (x)  + Vt12  (x)  t21  (x)  0 


P (x)T  (x)  P (x)  = 


(X)  Vt12  (x)t21  (x) 


on  0 < x < 1. 

Next  we  compute  G(x)  for  -l<x<0.  Since  G(x)  = G(-x) 
for  -1 < x < 0,  we  find  that 

G(x)  = [det  (A  (-x)  ) ]-1T  (-x)  -l<x<0. 

Moreover,  since  or  (x) , g (x) , y (x)  and  6 (x)  are  real,  the  rows 
of  P (x)  are  also  independent  eigenvectors  of  T(x)  for 
0<x<l.  Hence  P (-x)T  (-x)  P (-x)  1 is  diagonal  for  -l<x<0, 

A 

and  if  P(x)  is  defined  by 


P(x)  0 < x < 1 P(x)  0 < x < 1 

P (x)  = or  . 

P(-x)  -1  < x < 0 P(-x)  -1  < x < 0 


then  P(x)G(x)P~^  (x)  is  diagonal  for  -l<x<l.  if  in  addition 
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P(x)  can  be  extended  to  a function  P(z)  meroraorphic  in 
with  a pole  at  infinity,  then  E(z)  =P(z)  * (z)  defines  a 
sectionally  analytic  function  E(z),  satisfying  properties 
analogous  to  3.4  and  the  boundary  condition 


E (x)  = -D(x)  E (x)  +k(x) 


(3.19) 


where 


D(x)  * P(x)G(x)P_1(x) 


and  k(x)  = P(x)g(x).  An  example  illustrating  this  analysis 
appears  in  the  fourth  section. 


§4  Examples . 


As  a first  example,  we  consider  the  set  of  equations  (2.1), 
(2.2)  with  a(x)=a,  p(x)=p,  \(x)  = y and  6 (x)  = 6 . Here,  or, 
p,  Y»  and  6 are  all  complex  constants.  Since  this  example 


appears  in  various  applications,  we  display  the  appropriate 
matrices  and  write  out  the  solution.  We  assume  that 


The  coefficient  matrices  needed  in  (2.13)  and  (2.14)  are 


given  by: 


A-1  (x) 


1 

*6e2MTTi- 


(4.1) 


G(x)  = 


2(irri 


0f6  - py  - 2ig6  sin  nrr 
-Y0  L -2iay  sin  prr  c*6  - gy 


(4.2) 


P (x)  = -2i  sin  (an 


60 

V»0  Y6 

6p 

--v/apY6 

(4.3) 


while  the  diagonal  matrix,  D(x),  used  in  (2.15)  takes  the  form. 


D(x)  = 


«,6e2Kni-Ve 


<*6 -p y + 2i  sin  urrJatQyb  0 

0 a5  - p y - 2i  sin  iarr«/ttp 


(4.4) 


The  problem  is  then  reduced  to  the  solution  of  the  uncoupled 
Riemann  Hilbert  problem:  determine  Z^(z) , and  E2  (z)  analytic 
everywhere  in  the  plane  except  along  the  cut  (a,b)  where 
and  E2  satisfy  the  conditions. 


<(*)  = -e(vnU  f z;<x)  +V*)  . 


(4.5) 


E*(x)  = 
2 


- e(X'v)i  * zT  (x)  +k  (x)  , 

02  2 


(4.6) 


E1(z)  = E2(z)  = 0(|ZrA)  as  z- 


and  v*l«P*o»k^(x)  and  k2(x)  are  given  by: 
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v = tan 


* 

1 2 s in  inWa'g  v6 

. 6 or  - g v 


(4.7) 


X = tan-1  - tan 

2 2 *5 

p = [ (6 or  - 0y)  + 4 sin  utt(<*0Y6)] 


(4.8) 


(4.9) 


2 2 H 

v = [(O'6)  + (yp)  - 2of6  yp  cos2(jtt] 


(4.10) 


where 


K (x)  = 


kx  (x) 


k2  (x) 


[p(x)lrtx)] 


Fx(x) 


F2  <*> 


(4.11) 


For  example,  the  solution  of  the  coupled  pair  of  equations 


px  cp^tjdt  pl  cp2(t)dt 

J (x  - t)  + J (t  - x)  ~ fl  ' 


x € (-1,1) 


(4.12) 


pl  (t)  dt  px  cp2  (t)dt 

Jx  (t  - x/5  + Jq  (t  -x)^=  f2 


x € (-L  1) 


(4.13) 


(with  f^  and  f2  constant) 

require  the  determination  of  the  sectionally  analytic  functions 
E-^z),  I2(z)  with  E^  and  E2  vanishing  at  ® and  satisfying 
the  boundary  conditions. 


Ex(x)  = - E-^x)  +kx. 


E2(x)  = E2(x)  +k2. 


xe(-U) 


X €(-1,1) 


(4.14) 


(4.15) 


where  k^  and  k2  are  the  complex  constants 


I "1 
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kj=  (1  - i)  (fx+  f2),  k2=  - (l  + iM^-f^  . 


It  is  a simple  matter  to  demonstrate  that  the  solutions  to  the 
above  Riemann  Hilbert  problems  are  given  by 

El(z)  = 2nl(1  + l)  (fl+  V 0(z  +1)  (z-  1)  " 1}  (4,16) 


and 


E2<*>  - - 2k'1  + » «1-  f2>  ln  [iTl] 


(4.17) 


so  that 

$1  (z)  = - i[ (z)  + £2  (z)  ] 

and 

*2  (z)  = -i[£2(z)  - Ex(z)]  . 

The  unknown  functions  qp^  and  cp2  are  then  determined  from 
the  Abel  equations  (2.7)  and  (2.8)  where  R(x)  = 4</(l  - x^) . 

One  can  see  that  even  in  the  simple  case  (4.12)  and  (4.13)  the 
results  are  quite  complicated.  Fortunately,  for  applications, 
one  is  usually  only  interested  in  ^(z)  and  $2(z).  For 
problems  in  the  theory  of  elasticity,  [see  Lowengrub  [J  ]]  the 
stress  field  is  expressed  in  terms  of  these  two  sections lly 
analytic  functions. 

Secondly,  we  consider  the  particular  case  of  (3.1)  and 
(3.2)  with 

U = *S.  o(x)=p,  g(x)  = ^,  Y(x)=o-  and  6 (x)  = 


where  a and  g are  constants.  We  show  that  our  methods 
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yield  the  same  result  as  in  Lowengrub  [3],  We  consider  the 
coupled  system 

A(x)$+(x)  = - B(x)$“(x)  + F(x),  -1  < x < 1 


where  A (x)  , B (x) , and  F(x)  are  given  in  (3.14a)  while  A,  B 
* 

and  F are  defined  by  (3.16).  It  follows  that  the  condition 
det  A(x)/0  implies  that  p1 2-  a2 t 0.  In  addition,  if  0<x<l 
we  have 


G(x)  = - 


2 2 
0 - or 


2.  2 


X / X 


2iorB 


-2iag 


(^)J 


(4.18) 


-1 


P(x)  = 


(4.19) 


1 I (P~Q')“  0 

P (x)T  (x)P  (x)  = - - 

X 


(P+Qf)' 


while  if  -1  < x < 0 , 


(4.20) 


(p  - a) 


P(-x)  T(-x)  p-l(_x)  x 


(p  + a) 


2 (4.21) 


I,  2 2, 


and  det  A(-x)  = J(p  - a). 


Thus,  the  uncoupled  Riemann  boundary  value  problem  becomes: 

determine  two  sectionally  analytic  functions  £,  (z)  and  Z- (z) 

••  * 


that  vanish  at  infinity  and  along  the  cut  -1  ^ x £ 1 satisfy 
the  boundary  conditions 


zi(x)  = (p"+o)  rii(x)  + kl(x)  ' -1  < x < 1 

(4.22) 

E2(X)  = (pT?)  Vx)  +k2(x)  ' -Kx<l 

This  is  in  complete  agreement  with  Lowengrub  [.?  ],  where  the 
solution  to  the  original  pair  of  Abel  equations  is  given  along 
with  an  application  to  elasticity. 

§5-  Application  to  dual  relations. 

In  what  follows,  use  will  be  made  of  certain  well  known 
operators  of  fractional  integration,  differentiation  and  Hankel 
transforms  which  were  introduced  by  Erdelyi  and  Sneddon  [(  ]. 

Let  S denote  the  Hankel  transform 

T\f  Of 

co 

= 2“x"“Jo  ?1‘°'t<5>J'2t)+1,<S*>ai  . (5.D 

and  for  or  > 0 define  the  fractional  integral  operators  I 

r].a 
by 

, -2«*-2nfx  ~ ot-1  . 

-TO  V*-5’  s”  (5-2> 


f (?)d|. 


(5.3) 


K 


,„(£(!)  =x}  -££_  J <?2-xV  "5 


2ti, 


2 . a 17-2a-2n+l 


For  a < 0 the  fractional  derivatives  I and  K are 

T1»a  Tl»a 

defined  to  be  the  formal  inverses  of  the  operators  I 

ri+a,  -or 

and  KT]+a/ _a  respectively.  Frequent  use  will  be  made  of  the 
well  known  identities  relating  these  operators.  [See  Sneddon 
[ 7 ] , page  274j 

We  shall  consider  systems  of  dual  relations  of  the  form 


J [aA(  §)  +bB(?)]§“2QJ(j(5x)d§  = f x (x ) 0 < x < 1 

0 

00 

r [CA(?)  +dB(§)]§_2ej  (§x)d§  = f (x) 

J0  v 1 

00 

r a<s)J  <5x)a§  - o 
o p 

1 < X < ® 

B(s)J  (gx)ds  = 0 , 


<x> 


(5.4) 


where  a,b,c  and  d are  constants. 

Systems  similar  to  (1)  have  been  discussed  by  various 
authors.  Closed  forms  solutions  have  been  obtained  by  Lowengrub 
and  Sneddon  [ f ] for  the  cases  a = p = |i  = 0,  v = 1 and 

a = p = -h,  = v = ~H*  In  the  former  instance,  the  system 

(1)  was  reduced  to  a system  of  Carlemann  type  singular  integral 
equations  whereas  in  the  latter  the  generalized  Abel  system 


(3.1,  3.2)  was  obtained.  Westmann  used  similar  techniques  to 
construct  solutions  for  a = 0 and  n = v = 2,  whereas  Erdogan 
[ 2 ] presented  a method  for  reducing  (1)  to  an  infinite  system 
of  algebraic  equations. 

The  procedures  developed  for  treating  special  cases  of  (5.4) 
have  mostly  been  ad  hoc.  We  shall  indicate  a more  systematic 
approach  for  analyzing  (5.4)  which  in  certain  cases  reduces  to 
the  generalized  Abel  system  (3.1,  3.2). 

In  operator  notation  the  system  (5.4)  becomes 


J 


„ £ C aco  ( § ) +by(§)]?x3  =F.  (x), 
-a, 2a  1 


0 < x < 1 


(5.6) 


S { C cep  ( § ) + bf(?)],*x}  =F_(x),  0<x<l  (5.7) 

Y-S.2B 


S {cp(?);x}  = g.  (x)H  (1  - x) , x>l  (5.8) 

*0 


{*  (?) ;x} 


= g2  (x)H(l  - x)  , 


x > 1 


(5.9) 


2 2 ot 

where  cp(?)=A(?)/§  , Y(?)=B(?)/5  , F1(x)=2  "f^xj/x  , 

F2 (x)  = 22pf2 (x)/x2p  and  g1 (x)  and  g2  (x)  are  unknown 
functions  for  0<x<l. 


Formal  inversion  of  (5.8)  and  (5.9)  yields 
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= s {g,(x)  ? 5) 

u o 

= S oK  [k.  ] 

$<° 

- 2vrvJ  *l'vjM+Y(xt)k1  (x). 


(5.10) 


Y(§)  = S {g„ (x) ; ?} 

fr°  2 


s »K  [k  ] 

X0  £,6  2 


S JV 

^,6  2 

2 


6 -6  r 1-6  -r- 

= 2?  \x  k2(x)^+6(x?)dx 


(5.11) 


k.  (x) 

= K 

[g.  ] 

1 

^Y, 

*-Y  1 

k (x) 

= K 

Eg,] 

2 

£+6, 

-6  2 

and  v and  6 are  parameters  to  be  specified  later.  The 
manipulations  involved  in  (5.10)  and  (5.11)  are  well  known  and 
may  be  found  in  Sneddon  [ 7 J • 

Substitution  of  (5.10)  and  (5.11)  into  (5.6)  and  (5.7), 
followed  by  an  application  of  appropriate  fractional  transforms, 
yields 
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I S {Cacp(?)  +bY(?)];x]  = F (x) 

+ or,  \ - of,  2q ' A 


= a S CS  [k-]+bS  .S  [k  ] (5.12) 

^-a,X+2a  U,y  1 ^-afX+2a  X6 

fa  m £*  C* 


I ®s  [[ccp(?)  +dy(5)];x}  = F0  (x) 

*+&P  2 


= C S 


S [k,  ] +d  S o S [k_] 


g/P+2p  ^r.Y  ^-p»2p  ^.,6 


where  F (x)  = I [F  ], 

2+a'X 


F„(x)  = I [F,] 

2 ^ + p,  P 2 

2 


(5.13) 


and  X and  p are  parameters  as  yet  unspecified. 

There  are  two  trivial  cases  of  (5.12)  and  (5.13). 

If  |i  = v»  let  V = 6 , X = -o/  and  p=-p,  then  these  relations 


become 


S,  [k.  ] + b S„  . S„  [kj  = F.  (x) 


■'  ' w L Jfc-  J p fc/  W - W I J 4.  - 

^-a,a  ^.y  1 2~a,a  2'V  2 1 


= a K 


^ - a ,a+y 


[k  ] +b  K [k] 

i u , 2 

%-a,a+y 


(5.14) 


c S »S  [k_ ] + d S .S  [k_]  * 

■li_oo  -^v  ^ -U  - 0 0 -li  v 2 

2 P'P  2'y  2 P'P  2'Y 


F2(x) 


= C K [k. ] + d K Ck 

^ - b*b+v  1 rp-8+v 


(5.15) 
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If  we  invert  relations  (5.14)  and  (5.15),  one  obtains  a simple 
algebraic  system  for  k^  and  k2 

A second  trivial  case  results  if  + a - g = 0,  since  by 

choosing  X+o'  = p+  g = v + Q'  = 6+  e=  0 we  obtain 

a S °S  [kn]  + b S oS  [k_ ] = F.  (x) 

Ll  il  l U v21 

2 - a.  O'  2*  ~a  2 ~ a,a  2' 


= ak..  + b I [k„] 

L \i  £ 

2*a-P 

and 

c S S + d S S [k_]  = F_  (x) 

f-fS.B  2 * e,s  2'-» 

•cl  [k^+dkj. 


(5.16) 


(5.17) 


Application  of 


c k,  + d I 
1 


2' 


[k,] 

or-p 


to  (5.17)  yields 


(5.18) 


It  is  now  a simple  matter  to  solve  (5.16)  and  (5.18)  for  k^ 
and  k2. 

In  general  the  system  (5.12)  and  (5.13)  cannot  be  solved 
so  easily.  However,  it  may  be  simplified  in  either  of  two  ways. 
If,  for  y,  p and  6,  we  choose  y = -a,  6 = -g , p = - g 

and  X = - a then  the  system  becomes 
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a I [k.]  + b K [k  ] = F, 


(5.19) 


c K 


2 ~ 5*^+  (p-«) 


[k  ] +d  I [k  ] = F . 

1 \j  U— m 2 2 


U-v 
2'  2 


(5.20) 


Whereas  by  choosing  Y = 'k~^-0,,  5 = - a,  p = -g  and 


X = -a  we  obtain 

a K [k  ] + b I [k_]  = F, 

1 2 1 


(5.21) 


cl  [k.  ] + d K !>„]  = F„  . 

1 


(5.22) 


The  systems  (5.19)  - (5.22)  may  be  regarded  as  generalized 
Abel  systems.  In  most  applications  a = g and  hence  we  shall 
consider  only  this  case  for  the  remainder  of  this  section. 

It  should  be  observed  that  when  v - U is  an  even  integer  the 
operators  appearing  in  (5.19)  - (5.22)  are  not  of  fractional 
order.  It  is  then  a simple  matter  to  reduce  either  of  the 
systems  (5.19),  (5.20) 

or  (5.21),  (5.22)  to  a single  linear  ordinary  differential 

equation.  In  particular,  for  the  case  considered  by  Westmann  CSl 
(i.e.  v = |i  + 2)  the  system  (5.19),  (5.20)  reduces  to  a simple 

first  order  linear  differential  equation. 

If  (j  - v is  not  an  even  integer  then  the  operators  are 


_ - i_. 
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of  fractional  order,  with  both  fractional  integrals  and 
fractional  derivatives  appearing  in  both  systems.  However,  in 
certain  cases  it  is  possible  to  reduce  the  systems  (5.19)  - (5.22) 
to  Abel  systems  of  the  form  (3.1)  and  (3.2).  We  might  note 
that  if  p - v is  one  . the  Abel  system  obtained  will 

contain  only  operators  of  order  1/2,  whereas  in  general,  operators 
of  different  order  will  occur  within  the  same  system.  Provided 
neither  of  the  unknown  functions  appears  in  operators  of  differ- 
ent order,  a straight  forward  extension  of  the  technique 
presented  in  section  3 will  treat  such  systems. 

As  an  example,  consider  the  system  (1)  with  v = 1,  (i  = 0, 

<x  = (}  = h.  If  we  let  \ = p = 6 = -1  and  y = 0 we  obtain 


a K-%,/V+b  - *i<x) 


c a = f2(x) 


(5.23) 

(5.24) 


S ince , 


K0,-/*2^  "I  ^ ~x  * dt  Cfck2  (t)  3dt, 


*1 (t)  = g1(t)  and 


k2(t)  = K-is,lCg23  = 2t‘1ftg2(u)du. 


ml  x 


the  system  (5.23),  (5.24)  yields  the  generalized  Abel  system 


J (t2-*2)  i5tg1(t)dt-b  J (x2-t2)  ^g2  (t)  dt  =£^l)F1(x)  (5.25) 


0 < x < 1 


x 1 

X J (x2-fc2)  ijtg1(t)dt  + d J (t2-x2)_i5g2  (t)dt  = ^^(x)  (5.26 

0 x 

0 <x  < 1 . 
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